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Abstract 

Transferring quantum information between two qubits is a basic requirement for many applica- 
tions in quantum communication and quantum information processing. In the iterative quantum 
state transfer (IQST) proposed by D. Burgarth et al. [Phys. Rev. A 75, 062327 (2007)], this is 
achieved by a static spin chain and a sequence of gate operations applied only to the receiving 
end of the chain. The only requirement on the spin chain is that it transfers a finite part of the 
input amplitude to the end of the chain, where the gate operations accumulate the information. 
For an appropriate sequence of evolutions and gate operations, the fidelity of the transfer can 
asymptotically approach unity. We demonstrate the principle of operation of this transfer scheme 
by implementing it in a nuclear magnetic resonance quantum information processor. 
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I. INTRODUCTION 



Quantum state transfer (QST), i.e., the transfer of an arbitrary quantum state q;|0) + 
(3\1) from one qubit to another, is an important element in quantum computation and 
quantum communication . The most direct method to implement QST is based 

on SWAP operations This approach consist of a series of SWAP operations between 
neighboring qubits until the quantum state arrives at the target qubit. In a general-purpose 
quantum register, these quantum gates require the application of single- as well as two qubit 
operations. For longer distances, the number of such operations can become quite large; 
it may then be advantageous to rely on quantum teleportation instead j?], which requires 
fewer gate operations, but shared entanglement between sender and receiver. 

For specific systems, it is possible to transfer quantum information without applying gate 
operations, but instead relying on a static coupling network . The main difficulty with 
this approach is the required precision with which the couplings have to be realized in order 
to generate a transfer with high fidelity. 

This requirement can be relaxed significantly, without compromising the fidelity of the 
transfer, by applying gate operations to the receiving end of the spin chain that effects 
the transfer |5|. The capability for applying such gate operations is not an additional 
requirement, since such operations are required anyway if the spin chain is to be used 
for communication between quantum registers. This gate accumulates any amplitude of the 
initial state that is transferred along the chain. The protocol allows one, in principle, to 
obtain unit fidelity for the transfer, even if the couplings along the chain have arbitrary 
fiuctuations, as long as a finite amplitude reaches the end of the chain. Obtaining a large 
transfer amplitude requires multiple iterations, each of which includes the evolution of the 
spin chain and the two-qubit gate operation. The fidelity for transfer increases with the 
number of the iterations and can approach 1 asymptotically. Hence we refer to this protocol 
as the iterative quantum state transfer (IQST). In this paper we implement the protocol in 
an NMR quantum information processor and demonstrate its basic feasibility. 
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II. ITERATIVE TRANSFER ALGORITHM 



A. System 

We illustrate the IQST proposed in Ref. Q using a system of three spins coupled by 
Heisenberg XY- interactions, as shown in Figure [TJ The spin chain consists of spins 1 and 
2, which are coupled by a constant (time- independent) interaction. Spin 3 is the target spin 
used to receive the transferred quantum state. The interaction between spins 2 and 3 can 
be switched on and off. Our purpose is to transfer an arbitrary quantum state a\0) + f3\l) 
from spin 1 to 3, where a and (3 are two complex numbers normalized to|Q;p + |/5p = l. 

The Hamiltonian of the the spin chain without the end qubit is 

H^2 = lriJu{cTlal + ayy), (1) 
where Ju denotes the coupling strength. The Hamiltonian of spins 2 and 3 is 

H23{t) = lnJ23m^lcrl + alal), (2) 
where J23{t) is J23 when the interaction is switched on and otherwise. 



B. IQST algorithm 

The purpose of the IQST algorithm is the transfer of an arbitrary state a\0) + from 
the start of the chain (qubit 1) to the end (qubit 3). We start the discussion by choosing 
as the initial state of the complete 3-qubit system the state q;|000) + /?|100), i.e. a product 
state with spin 1 in state a|0) + and spins 2 and 3 in |0). Transferring the a|0) part of 
the input state is trivial, since spins 1 and 3 are in the same state and this state is invariant 
under the XY interaction. We therefore only have to consider the (3\1) part. 

The chosen initial state of the spin chain is not unique. We could, e.g., choose to start 
with the total system in a|011) + In this case, the |111) is invariant and only the 

transfer of the a|0) part needs to be considered. At the end of this section, we discuss 
additional possibilities. 

The iterative transfer scheme of Burgarth et al. consists of a continuous evolution under 
the spin-chain Hamiltonian, interrupted by successive applications of the end-gate operation. 
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We write the transfer operator as 



(3) 
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represents the evolution of the spin chain and 
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the end gate operation. Here, C12 = cos(7rJi2T) and 5'i2 = sin(7rJi2T) and n represents the 
iteration step. The parameters c^, d„ are related by the unitarity condition |c„p+ = 1. 
For each step of the iteration, they are equal to the coefficients of the relevant states |010) 
and 1 001) just before the gate is applied. Under this condition, 

W^%Cn,dn)icn\OlO) + dn\00l)) = |001), 

i.e. the transfer to the final state |001) is maximized. 
During the n*'* step, the two coefficients are 
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C. Quantification of transfer 



After k iterations, |100) is transferred to 



l^fe) = TfcllOO) = Cf^llOO) + Jl- CfllOOl) 



(8) 
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Apparently, the transfer increases monotonically with the number of iterations and can 
asymptotically approach unity provided IC12I < 1. Writing Fk = (001|^'a;) for the overlap of 
the system with the target state, we find 

Fk = Vl - C!l. (9) 

Eq. ([3]) implies that only the spin chain or the end gate are active at a given time. If the 
spin chain interactions are static (not switchable), this can only be realized approximately 
if the coupling between the two end-gate qubits is much stronger than the couplings in the 
spin chain, J23 ^ Ji2- In the NMR system, we instead refocus the spin-chain interaction 
during the application of the end-gate operation to better approximate the ideal operation 

where 

tan{7lJ23tn) = -iCn/dn- (H) 



D. Generalization to mixed states 



The IQST algorithm works also when the spin chain is in a suitable mixed state. As 
an example, we choose a = P = The second and third qubit can be chosen in any 
combination of |0) and Here, we implement all four possibilities in parallel by 
putting qubits 2 and 3 into the maximally mixed state ® J^, where I denotes the unit 
operator and the upper index labels the qubit. The sample thus contains an equal number 
of molecules with qubits in the states a\Ol) + 0\ll) with / = {00,01, 10, 11}. The traceless 
part of the corresponding density operator is 9j 

P^n^=J:CTl^i\l){l\). (12) 
1=00 

If the system is initially in one of the states \l) = |01), |10), it acquires an overall phase 
factor of —1 during the transfer. Combining this with the results of Sec. IIIBl we find that 
after k iterations, the system is in the state 



Pk = n Prm Tl = Jl- Fl al I' I' + Fk alalal (13) 
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Similarly, when the initial state is chosen as 

P^n^=J:crl(S{\l){l\), (14) 

the algorithm generates the state 



y 

1=00 



Pk = T,p,,mTl = ^jl- Fl all'l' + F, a\alal (15) 

after k iterations. 

III. IMPLEMENTATION 

For the experimental implementation, we chose the ^H, ^^F, and ^"^C spins of Ethyl 2- 
fiuoroacetoacetate as qubits. The chemical structure of Ethyl 2-fluoroacetoacetate is shown 
in Figure [21 where the three qubits are denoted as HI, F2, and C3, respectively. The 
strengths of the J-couplings are J\2 = 48.5 Hz, J23 = —195.1 Hz and J13 = 160.8 Hz. Ti 
and T2 values for these three nuclei are listed in the right table in Figure [2l In the rotating 



frame, the Hamiltonian of the three- qubit system is 



10 



111 



Unmr = 2 (-^i2^>' + J23(rlcTl + Jncrlal). (16) 

The sample consisted of a 3:1 mixture of unlabeled Ethyl 2-fluoroacetoacetate and d6- 
acetone. Molecules with a ^'^C nucleus at position 2, which we used as the quantum register, 
were therefore present at a concentration of about 1%. They were selected against the 
background of molecules with ^^C nuclei by measuring the ^^C signal. We chose HI as the 
input qubit and C3 as the target qubit. Figure [3] (a) shows the ^^C NMR spectrum obtained 
by applying a readout pulse to the system in its thermal equilibrium state. Each of the 
resonance lines is associated with a specific spin state of qubits 1 and 2. 



A. Initial state preparation 

The initial pseudo-pure state |000) is prepared by spatial averaging (l2l|. The following 
radio- frequency (rf) and magnetic field gradient pulse sequence transforms the system from 
the equilibrium state 

Peg = + 12(^1 + (17) 
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to 1000): - [HI - [gradl - [n/2]l - [1/2 J13] - [-vr/2]^ - [vr/4]^ - [-1/2J,;] - [-vr/4]^ - 
[grad]z — [7r/4];[. — [I/2J13] — [— 7r/4]^ — [grad]^. Here 71, 72 and 73 denote the gyromagnetic 
ratios of HI, F2, and C3, respectively, and cos0i = 273/71, and cos02 = 73/272- [grad]z 
denotes a gradient pulse along the z- axis, [vr/2];^ denotes a 7r/2 pulse along the x- axis 
acting on the HI qubit. Overall phase factors have been ignored. 

The coupled-spin evolution between two spins, for instance, [1/2 J13], can be realized by 
the pulse sequence l/4Ji3 — [- 7r]| — l/4Ji3 — [— Tr]^, where l/4Ji3 denotes the evolution caused 



by Hnmr for a time 1/4 J13 [13 1. 



The target state can be prepared directly from the state |000) by applying a [vr/2]^ pulse. 
It corresponds to |00)(|0) -\l))/^/2, i.e. to transverse magnetization of the target spin, with 
the first two qubits in state |00). If we measure the free induction decay (FID) of this state 
and calculate the Fourier transform of the signal, we obtain the spectrum shown in Figure 
[3] (b). This spectrum serves as the reference to which we scale the data from the IQST 
experiment. 

The input state for the IQST is \^in) = \'ip{6))\QQ) . We generate this state by rotating 
HI by an angle around the y-axis: = e^^'"i/2|ooo). After k iterations of the IQST 

algorithm, \^in) is transferred to 



Tfc|^,„) = [(1 - Ffc) cos(0/2)|O) -Jl- F2sin(^^/2)|1)]|00) + |00)F,#(^)). (18) 



Here, we have used Eqs. (|8l|9]) and assumed C12 > 0, without loss of generality. Hence the 
state transfer can be observed through measuring carbon spectra. 

For the mixed input state, [Eq. f|T^ ] can be generated from p^q through the pulse 

H 



sequence 



7r,3 r ,1 rTT 



W.-y- (19) 



B. Effective XY- interact ions 



The IQST algorithm requires XY interactions, while the natural Hamiltonian contains 
ZZ couplings. To convert the ZZ interactions into XY type, we decompose the evolution 
Q-ivi'^x'^i+'^l'^y) into e~*'^'^^'^^e~*'^°"S'^!' [l5|] using [c^cr^, 0"yCry] = 0, where ip denotes an arbitrary 
real number. These tranformations can be implemented by a combination of radio-frequency 
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pulses and free evolutions under the J-couplings: 



16| 



g-i(pcrJo-i ^ g±j7ro-*/4g±i7r4/4g-i(pcrJ4gq:i7ro-J/4g:fi7ro-J^/4 ^20) 
^-iipa^cj'y ^ g±j7ro-*/4g±i7ro-i;/4g-j</3crJ4g=Fi?r(T^/4g=Fi?i-o-i/4 _ ^21) 

Figure H] shows the complete pulse sequence for implementing the IQST, starting from 
The subscript n indicates that the pulses in the square brackets have to be repeated for 
every iteration. The duration of each W^^ segment varies, since tn = — arctan(ic„/c/„)/7r J23. 

For the initial state pmi in Eq. ( |T2l) . the propagators n can be simplified: since the 
density operator commutes with alal and cr^cr^ at all times, it is sufficient to generate the 
propagator 

g-j7rJ23tnO-xO"S/2g-i7!-Jl2-ro-Jo-2/2 



Similarly, for the initial state in Eq. flMj) . iteration n can be replaced by 
g-«7rj23tno-acry/2g-i7rJi2ro-J(T2/2_ |;]^ggg simplified verslous to shorten the duration of the 

experiment and thereby increase the fidelity. 



C. Results for state transfer 



When r = 1/2 J12, the transfer can be implemented in a single step with a theoretical 
fidelity of 100%. The state transfer from HI to C3 can be observed by measuring ^^C spectra. 
The experimental result for |\l/i„) = |?/'(7r/4))|00) is shown in Figure [5] (a). Comparing with 
Figure [3] (b) one finds that the output state is |00)(|0) — |l))/-\/2, i.e., the state \ip{7r/A)) is 
transferred from HI to C3. 

Figure [5] (b), show the corresponding result for the transfer of ay from HI to C3 in a single 
step, with qubits 2 and 3 initially in the completely mixed state. For this experiment, the 
receiver phase was shifted by 7r/2 with respect to the upper spectrum. Since this experiment 
implements the transfer for all possible states of the other qubits in parallel, we observe four 
resonance lines corresponding to the states {00, 01, 10, 11} of qubits 1 and 2. For the states 
with odd parity, the transfer adds an overall phase factor of -1, which is directly visible as 
a negative amplitude in the spectrum. 

To demonstrate that iterative transfer works for a range of coupling strengths or (equiva- 
lently) evolution periods, we chose r = 1/5 J12 and r = 1/6 Ju- For the case of pseudo-pure 
input states, three iterations are implemented for either case. When 6 changes from to 27r 
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the experimental results obtained from these transfer experiments are summarized in Figure 
[6l where the vertical axis denotes the amplitude of the NMR spectrum. For each input 
state the amplitude increases with the number of iterations. The increase of the amplitude 
shows the increase of the fidelity for the state transfer. The dependence on the input state 
parameter 9 has the expected sin(^) dependence. 

The experimental data obtained for the mixed input states are summarized in Figures [7] 
(a) and (b), for r = 1/5 J12 and r = 1/6 J12, respectively. The positive lines indicate that the 
transfer occurs with positive sign if qubits 1 and 2 are in state |00) or |11), and with negative 
sign for the states |01) or |10), in agreement with Eq. ( |T5l) . Obviously the amplitude of the 
signals increases with the number of iterations. According to Eq. fllSp the increase of the 
amplitudes is a direct measure for the progress of the quantum state transfer. 

IV. DISCUSSION AND CONCLUSION 

Our results clearly demonstrate the validity of the iterative state transfer algorithm of 
Burgarth et al. In principle, it is possible to iterate the procedure indefinitely, always 
improving the fidelity of the transfer. In practice, every iteration also increases the amount 
of signal loss, either through decoherence or through experimental imperfections. 

According to Eq. f[T^ . the fidelity of the transfer is 

F, = \Tr[{alalal)pk]\. (22) 

The experimental measurement corresponds to a summation of the amplitudes of the reso- 
nance lines. We normalized the experimental values to the amplitudes obtained by direct 
preparation of the target states [see Figure [3] (a)]. In Figure [H we show the experimen- 
tally measured fidelities of the transfer of the state ay for 1-5 iterations. As expected, the 
experimental data points are below the theoretical curves (full lines). 

The experimental points can be fitted quite well if we include a decay parameter for 
each iteration. The dashed curves in Figure [H] represent the function Fke~^^ with r = 0.087 
and r = 0.079 for r = 1/5 J12 and r = 1/6 J12, respectively. Each iteration thus adds 
imperfections (experimental plus decoherence) of about 8 %. Larger numbers of iterations 
are meaningful only if this error rate can be reduced. 

In conclusion, we have implemented the iterative quantum state transfer in a three qubit 
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NMR quantum information processor. The result shows that it is indeed possible to accu- 
mulate the quantum state at the end of a Heisenberg spin chain, whose couplings are always 
active. 
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FIG. 1: The spin chain including the target spin (3) used for implementing the IQST. The XY- 
interactions in the spin chain, denoted by the solid line, is always active, while the XY- interaction 
between spins 2 and 3, denoted by the dashed line, can be switched on and off. W'^^ denotes the 
end gate applied to spins 2 and 3. U^'^ denotes the evolution of spin chain. 
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FIG. 2: (Color online) The chemical structure of Ethyl 2-fluoroacetoacetate. The three spins in the 
dashed oval are the three qubits for implementing IQST. The strengths (in Hz) of the J-couplings 
between the relevant nuclear spins and the relaxation times are listed in the left and right tables, 
respectively. 
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FIG. 3: (a) -"^^C NMR spectrum obtained by applying a selective readout pulse to the system in its 
thermal equilibrium state. The four resonance lines correspond to specific states of the spin-chain 
qubits HI and F2, as indicated by the labels above the resonance lines. The assignment takes into 
account that J13 > and J23 < 0. (b) ^^C NMR spectrum of the state |00)(|0) - |1))/V2, which 
was obtained by applying a [7r/2]^ pulse to |000). 
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FIG. 4: (Color online) Pulse sequence for implementing the IQST. The two blocks that implement 
U^'^{t) and W'^^{cn,dn) are separated by the dash-dotted line and "[...]n" indicates iteration n. 
The delays tn are given by Eq. pT]) . The narrow rectangles denote 7r/2 pulses, and the wide ones 
denote vr pulses, where x, —x, y, or —y denote the direction along which the pulse is applied. The vr 
pulses are applied in pairs with opposite phases to reduce experimental errors p^]. The durations 
of the pulses are so short that they can be ignored. 
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FIG. 5: Experimental results for quantum state transfer with r = 1/2 Ju- The initial states are 
[|0) — |l)]|00)/\/2 and cr^, corresponding to figures (a) and (b), respectively. In the first experiment, 
the receiver phase was set to x, in the second experiment it was set to y. 
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FIG. 6: (Color online) Experimental results for demonstrating the IQST when the initial state is 
[cos(0/2)|O) — sin(0/2)|l)]|OO). Two cases for r = l/5Ji2 and r = I/6J12 are shown in Figures 
(a) and (b). For each case three iterations are implemented. The experimental data after the 
completion of iteration 1, 2, and 3 are marked by "*", "+", and "x", respectively. The data can 
be fitted as sin functions of which amplitudes represent the measured fidelities experimentally. The 
dashed curves show sin(0). 
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FIG. 7: (Color online) ^^C NMR spectra demonstrating the IQST of the state for r = 1/5 J12 
and r = 1/6 Ji2. For each case, the spectra after the completion of iteration 1, 2, and 3 are shown 
as the blue, black and red cm'ves, respectively. The resonance lines corresponding to the |00) state 
of the spin chain are enlarged in the inset. The dashed curves are the corresponding sections of 
the reference spectrum in Figure [3] (a). 
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FIG. 8: (Color online) Experimentally measured fidelity of the iterative state transfer as a function 
of the number of iteration steps when r = 1/5 J12 and r = 1/6 J12. The experimental data are 
fitted to the function FkC'^'' with r = 0.087 and 0.079 for the two cases, respectively. The two 
solid curves represent the theoretical fidelities for ideal conditions, and the two dash-dotted 
curves correspond to e~^^ . The dark and light curves correspond to the cases of r = l/5Ji2 and 
r = 1/6 J12, respectively. 



18 



